Gaussian accelerated molecular dynamics simulations facilitate prediction of the permeability of cyclic peptides.

Despite their widespread use as therapeutics, clinical development of small molecule drugs remains challenging. Among the many parameters that undergo optimization during the drug development process, increasing passive cell permeability (i.e., log(P)) can have some of the largest impact on potency. Cyclic peptides (CPs) have emerged as a viable alternative to small molecules, as they retain many of the advantages of small molecules (oral availability, target specificity) while being highly effective at traversing the plasma membrane. However, the relationship between the dominant conformations that typify CPs in an aqueous versus a membrane environment and cell permeability remain poorly characterized. In this study, we have used Gaussian accelerated molecular dynamics (GaMD) simulations to characterize the effect of solvent on the free energy landscape of lariat peptides, a subset of CPs that have recently shown potential for drug development (Kelly et al., JACS 2021). Differences in the free energy of lariat peptides as a function of solvent can be used to predict permeability of these molecules, and our results show that permeability is most greatly influenced by N-methylation and exposure to solvent. Our approach lays the groundwork for using GaMD as a way to virtually screen large libraries of CPs and drive forward development of CP-based therapeutics.

Lipid membrane mimetics and oligomerization tune functional properties of proteorhodopsin.

The functional properties of proteorhodopsin (PR) have been found to be strongly modulated by oligomeric distributions and lipid membrane mimetics. This study aims to distinguish and explain their effects by investigating how oligomer formation impacts PR's function of proton transport in lipid-based membrane mimetic environments. We find that PR forms stable hexamers and pentamers in both E. coli membranes and synthetic liposomes. Compared with the monomers, the photocycle kinetics of PR oligomers is ∼2 and ∼4.5 times slower for transitions between the K and M and the M and N photointermediates, respectively, indicating that oligomerization significantly slows PR's rate of proton transport in liposomes. In contrast, the apparent pKa of the key proton acceptor residue D97 (pKaD97) of liposome-embedded PR persists at 6.2-6.6, regardless of cross-protomer modulation of D97, suggesting that the liposome environment helps maintain PR's functional activity at neutral pH. By comparison, when extracted directly from E. coli membranes into styrene-maleic acid lipid particles, the pKaD97 of monomer-enriched E50Q PR drastically increases to 8.9, implying that there is a very low active PR population at neutral pH to engage in PR's photocycle. These findings demonstrate that oligomerization impacts PR's photocycle kinetics, while lipid-based membrane mimetics strongly affect PR's active population via different mechanisms.

Physicochemical Properties Altered by the Tail Group of Lipid Membranes Influence Huntingtin Aggregation and Lipid Binding.

Huntington's disease is a neurodegenerative disorder caused by an expanded polyglutamine (polyQ) domain within the huntingtin protein (htt) that initiates toxic protein aggregation. Htt directly interacts with membranes, influencing aggregation and spurring membrane abnormalities. These interactions are facilitated by the 17 N-terminal residues (Nt17) that form an amphipathic α-helix implicated in both lipid binding and aggregation. Here, the impact of unsaturation in phospholipid tails on htt-lipid interaction and htt aggregation was determined. There was no correlation between the degree of htt-lipid complexation and the degree of htt aggregation in the presence of each lipid system, indicating that lipid systems with different properties uniquely alter the membrane-mediated aggregation mechanisms. Also, the association between Nt17 and membrane surfaces is determined by complementarity between hydrophobic residues and membrane defects and how easily the peptide can partition into the bilayer. Our results provide critical insights into how membrane physical properties influence downstream htt aggregation.

Homo-oligomerization of the human adenosine A2A receptor is driven by the intrinsically disordered C-terminus.

G protein-coupled receptors (GPCRs) have long been shown to exist as oligomers with functional properties distinct from those of the monomeric counterparts, but the driving factors of oligomerization remain relatively unexplored. Herein, we focus on the human adenosine A2A receptor (A2AR), a model GPCR that forms oligomers both in vitro and in vivo. Combining experimental and computational approaches, we discover that the intrinsically disordered C-terminus of A2AR drives receptor homo-oligomerization. The formation of A2AR oligomers declines progressively with the shortening of the C-terminus. Multiple interaction types are responsible for A2AR oligomerization, including disulfide linkages, hydrogen bonds, electrostatic interactions, and hydrophobic interactions. These interactions are enhanced by depletion interactions, giving rise to a tunable network of bonds that allow A2AR oligomers to adopt multiple interfaces. This study uncovers the disordered C-terminus as a prominent driving factor for the oligomerization of a GPCR, offering important insight into the effect of C-terminus modification on receptor oligomerization of A2AR and other GPCRs reconstituted in vitro for biophysical studies.

Intramolecular interactions play key role in stabilization of pHLIP at acidic conditions.

The pH-Low Insertion Peptide (pHLIP) is a membrane-active peptide that spontaneously folds into a transmembrane α-helix upon acidification. This activity enables pHLIP to potentially act as a vector for drugs related to diseases characterized by acidosis such as cancer or heart ischemia. Presently, due to aggregation-based effects, formulations of pHLIP are only viable at near-µM concentrations. In addition, since most of pHLIP's measurable qualities involve a membrane, probing the details of pHLIP in the interstitial region is difficult. In attempts to shed light on these issues, we performed constant pH molecular dynamics simulations on pHLIP as well as P20G, a variant with increased helicity, in solution at 0 and 150 mM NaCl over a broad range of pHs. In general, the addition of ions reduced the effective pKa of the acidic residues in pHLIP. P20G exhibits a higher helicity than pHLIP in general and is more compact than pHLIP at pH values under 4. In terms of charge effects, sodium cations localized predominantly to the C-terminus of the peptide with a high density of acidic residues. Additionally, the salt bridge between R11 and D14 is by far the most favored and particularly so with pHLIP at 150 mM NaCl. We expect that this approach will be a valuable tool to screen variants of pHLIP for favorable properties in solution, an aspect of pHLIP design that to this point has largely been neglected.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides.

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Bacterial Analogs to Cholesterol Affect Dimerization of Proteorhodopsin and Modulates Preferred Dimer Interface.

Hopanoids, the bacterial analogues of sterols, are ubiquitous in bacteria and play a significant role in organismal survival under stressful environments. Unlike sterols, hopanoids have a high degree of variation in the size and chemical nature of the substituent attached to the ring moiety, leading to different effects on the structure and dynamics of biological membranes. While it is understood that hopanoids can indirectly tune membrane physical properties, little is known on the role that hopanoids may play in affecting the organization and behavior of bacterial membrane proteins. In this work we used coarse-grained molecular dynamics simulations to characterize the effects of two hopanoids, diploptene (DPT) and bacteriohopanetetrol (BHT), on the oligomerization of proteorhodopsin (PR) in a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phophoethanolamine (POPE) and 1-palmitoyl-2-oleoyl-sn-3-phosphoglycerol (POPG). PR is a bacterial membrane protein that functions as a light-activated proton pump. We chose PR based on its ability to adopt a distribution of oligomeric states in different membrane environments. Furthermore, the efficiency of proton pumping in PR is intimately linked to its organization into oligomers. Our results reveal that both BHT and DPT indirectly affect dimerization by tuning membrane properties in a fashion that is concentration-dependent. Variation in their interaction with PR in the membrane-embedded and the cytoplasmic regions leads to distinctly different effects on the plasticity of the dimer interface. BHT has the ability to intercalate between monomers in the dimeric interface, whereas DPT shifts dimerization interactions via packing of the interleaflet region of the membrane. Our results show a direct relationship between hopanoid structure and lateral organization of PR, providing a first glimpse at how these bacterial analogues to eukaryotic sterols produce very similar biophysical effects within the cell membrane.

Activation of the G-Protein-Coupled Receptor Rhodopsin by Water.

Visual rhodopsin is an important archetype for G-protein-coupled receptors, which are membrane proteins implicated in cellular signal transduction. Herein, we show experimentally that approximately 80 water molecules flood rhodopsin upon light absorption to form a solvent-swollen active state. An influx of mobile water is necessary for activating the photoreceptor, and this finding is supported by molecular dynamics (MD) simulations. Combined force-based measurements involving osmotic and hydrostatic pressure indicate the expansion occurs by changes in cavity volumes, together with greater hydration in the active metarhodopsin-II state. Moreover, we discovered that binding and release of the C-terminal helix of transducin is coupled to hydration changes as may occur in visual signal amplification. Hydration-dehydration explains signaling by a dynamic allosteric mechanism, in which the soft membrane matter (lipids and water) has a pivotal role in the catalytic G-protein cycle.

Using Simulation to Understand the Role of Titration on the Stability of a Peptide-Lipid Bilayer Complex.

The pH-low insertion peptide (pHLIP) is an anionic membrane-active peptide with promising potential for applications in imaging of cancer tumors and targeted delivery of chemotherapeutics. The key advantage of pHLIP lies in its acid sensitivity: in acidic cellular environments, pHLIP can insert unidirectionally into the plasma membrane. Partitioning-folding coupling is triggered by titration of the acidic residues in pHLIP, transforming pHLIP from a hydrophilic to a hydrophobic peptide. Despite this knowledge, the reverse pathway that leads to exit of the peptide from the plasma membrane is poorly understood. Our hypothesis is that sequential deprotonation of pHLIP is a prerequisite for exit of the peptide from the plasma membrane. We carried out molecular dynamics (MD) simulations to characterize the effect that deprotonation of the acidic residues of pHLIP has on the stability of the peptide when inserted into a model lipid bilayer of 1-palmitoyl-2-oleoyl-sn-3-phosphocholine (POPC). Initiation of the exit mechanism is facilitated by a complex relationship between the peptide, bulk solvent, and the membrane environment. As the N-terminal acidic residues of pHLIP are deprotonated, localized loss of helicity drives unfolding of the peptide and more pronounced interactions with the bilayer at the lipid-water interface. Deprotonation of the C-terminal acidic residues (D25, D31, D33, and E34) leads to further loss of secondary structure distal from the C-terminus, as well as formation of a water channel that stabilizes the orientation of pHLIP parallel to the membrane normal. Together, these results help explain how stabilization of intermediates between the surface-bound and inserted states of pHLIP occur and provide insights into rational design of pHLIP variants with modified abilities of insertion.

Theoretical Insights into the Mechanism of Wavelength Regulation in Blue-Absorbing Proteorhodopsin.

Proteorhodopsin (PR) is a light-driven proton pump that is most notable for ushering in the discovery of an ever-increasing number of microbial retinal proteins that are at the forefront of fields such as optogenetics. Two variants, blue (BPR) and green (GPR) proteorhodopsin, have evolved to harvest light at different depths of the ocean. The color-tuning mechanism in PR is controlled by a single residue at position 105: in BPR it is a glutamine, whereas in GPR it is a leucine. Although the majority of studies on the spectral tuning mechanism in PR have focused on GPR, detailed understanding of the electronic environment responsible for spectral tuning in BPR is lacking. In this work, several BPR models were investigated using quantum mechanics/molecular mechanics (QM/MM) calculations to obtain fundamental insights into the color tuning mechanism of BPR. We find that the molecular mechanism of spectral tuning in BPR depends on two geometric parameters, the bond length alternation and the torsion angle deviation of the all-trans-retinyl chromophore. Both parameters are influenced by the strength of the hydrogen-bonded networks in the chromophore-binding pocket, which shows how BPR is different from other microbial rhodopsins.

Ions Modulate Key Interactions between pHLIP and Lipid Membranes.

The pH-low insertion peptide (pHLIP) is used for targeted delivery of drug cargoes to acidic tissues such as tumors. The extracellular acidosis found in solid tumors triggers pHLIP to transition from a membrane-adsorbed state to fold into a transmembrane α-helix. Different factors influence the acidity required for pHLIP to insert into lipid membranes. One of them is the lipid headgroup composition, which defines the electrostatic profile of the membrane. However, the molecular interactions that drive the adsorption of pHLIP to the bilayer surface are poorly understood. In this study, we combine biophysical experiments and all-atom molecular dynamics simulations to understand the role played by electrostatics in the interaction between pHLIP and a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer. We observed that the solution ionic strength affects the structure of pHLIP at the membrane surface as well as the acidity needed for different steps in the membrane insertion process. In particular, our simulations revealed that an increase in ionic strength affected both pHLIP and the bilayer; the coordination of sodium ions with the C-terminus of pHLIP led to localized changes in helicity, whereas the coordination of sodium ions with the phosphate moiety of the phosphocholine headgroups had a condensing effect on our model bilayer. These results are relevant to our understanding of environmental influences on the ability of pHLIP to adsorb to the cell membrane and are useful in our fundamental understanding of the absorption of pH-responsive peptides and cell-penetrating peptides.

Cooperative Nonbonded Forces Control Membrane Binding of the pH-Low Insertion Peptide pHLIP.

Peptides with the ability to bind and insert into the cell membrane have immense potential in biomedical applications. pH (low) insertion peptide (pHLIP), a water-soluble polypeptide derived from helix C of bacteriorhodopsin, can insert into a membrane at acidic pH to form a stable transmembrane α-helix. The insertion process takes place in three stages: pHLIP is unstructured and soluble in water at neutral pH (state I), unstructured and bound to the surface of a membrane at neutral pH (state II), and inserted into the membrane as an α-helix at low pH (state III). Using molecular dynamics simulations, we have modeled state II of pHLIP and a fast-folding variant of pHLIP, in which each peptide is bound to a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer surface. Our results provide strong support for recently published spectroscopic studies, namely that pHLIP preferentially binds to the bilayer surface as a function of location of anionic amino acids and that backbone dehydration occurs upon binding. Unexpectedly, we also observed several instances of segments of pHLIP folding into a stable helical turn. Our results provide a molecular level of detail that is essential to providing new insights into pHLIP function and to facilitate design of variants with improved membrane-active capabilities.

Allosteric Effects of the Proton Donor on the Microbial Proton Pump Proteorhodopsin.

Proteorhodopsin (PR) is a microbial proton pump that is ubiquitous in marine environments and may play an important role in the oceanic carbon cycle. Photoisomerization of the retinal chromophore in PR leads to a series of proton transfers between specific acidic amino acid residues and the Schiff base of retinal, culminating in a proton motive force to facilitate ATP synthesis. The proton donor in a similar retinal protein, bacteriorhodopsin, acts as a latch to allow the influx of bulk water. However, it is unclear if the proton donor in PR, E108, utilizes the same latch mechanism to become internally hydrated. Here, we used molecular dynamics simulations to model the changes in internal hydration of the blue variant of PR during photoactivation with the proton donor in protonated and deprotonated states. We find that there is a stark contrast in the levels of internal hydration of the cytoplasmic half of PR based on the protonation state of E108. Instead of a latch mechanism, deprotonation of E108 acts as a gate, taking advantage of a nearby polar residue (S61) to promote the formation of a stable water wire from bulk cytoplasm to the retinal-binding pocket over hundreds of nanoseconds. No large-scale conformational changes occur in PR over the microsecond timescale. This subtle yet clear difference in the effect of deprotonation of the proton donor in PR may help explain why the photointermediates that involve the proton donor (i.e., M and N states) have timescales that are orders of magnitude different from the archaeal proton pump, bacteriorhodopsin. In general, our study highlights the importance of understanding how structural fluctuations lead to differences in the way that retinal proteins accomplish the same task.

Protonation Enhances the Inherent Helix-Forming Propensity of pHLIP.

Cell-penetrating peptides (CPPs) can be potentially used in targeted delivery of therapeutic cargoes. However, their conformation in solution is poorly understood. We employed molecular dynamics simulations to probe the structural fluctuations of an anionic CPP, pH Low Insertion Peptide (pHLIP), in solution to determine the effects of modifications to selected residues on the structure of pHLIP. Two types of modifications were tested: (1) protonation of aspartic acid residues and (2) point mutations known to affect the acid sensitivity of pHLIP. pHLIP samples conformations ranging from coil to helix to sheet, and modifications to pHLIP lead to subtle shifts in the balance between these conformations. In some instances, pHLIP is as likely to form a helical conformation as it is to form an unstructured coil. Understanding the behavior of pHLIP in solution is necessary for determining optimal conditions for administration of pHLIP and design of promising pHLIP variants.

Molecular Dynamics Simulations as a Tool for Accurate Determination of Surfactant Micelle Properties.

Molecular dynamics (MD) simulations were used to characterize the equilibrium size, shape, hydration, and self-assembly of dodecylphosphocholine (DPC) and dodecyl-ß-D-maltoside (DDM) micelles. We show that DPC molecules self-assemble to form micelles with sizes within the range reported in the experimental literature. The equilibrium shape of DPC and DDM micelles as well as associated micellar radii are in agreement with small-angle X-ray scattering (SAXS) experiments and theoretical packing parameters. In addition, we show that hydration of the micelle interior is limited; however, flexibility of the acyl chains leads to dynamic encounters with the solvated outer shell of the micelle, providing an explanation for long-standing differences in models of micelle hydration. Altogether, our results provide fundamental understanding of physical characteristics of micelles that can be utilized to study other types of detergents and proteomicelle complexes.

Structural Insights into Selective Ligand-Receptor Interactions Leading to Receptor Inactivation Utilizing Selective Melanocortin 3 Receptor Antagonists.

Systematic N-methylated derivatives of the melanocortin receptor ligand, SHU9119, lead to multiple binding and functional selectivity toward melanocortin receptors. However, the relationship between N-methylation-induced conformational changes in the peptide backbone and side chains and melanocortin receptor selectivity is still unknown. We conducted comprehensive conformational studies in solution of two selective antagonists of the third isoform of the melanocortin receptor (hMC3R), namely, Ac-Nle-c[Asp-NMe-His6-d-Nal(2')7-NMe-Arg8-Trp9-Lys]-NH2 (15) and Ac-Nle-c[Asp-His6-d-Nal(2')7-NMe-Arg8-NMe-Trp9-NMe-Lys]-NH2 (17). It is known that the pharmacophore (His6-DNal7-Arg8-Trp9) of the SHU-9119 peptides occupies a ß II-turn-like region with the turn centered about DNal7-Arg8. The analogues with hMC3R selectivity showed distinct differences in the spatial arrangement of the Trp9 side chains. In addition to our NMR studies, we also carried out molecular-level interaction studies of these two peptides at the homology model of hMC3R. Earlier chimeric human melanocortin 3 receptor studies revealed insights regarding the binding and functional sites of hMC3R selectivity. Upon docking of peptides 15 and 17 to the binding pocket of hMC3R, it was revealed that Arg8 and Trp9 side chains are involved in a majority of the interactions with the receptor. While Arg8 forms polar contacts with D154 and D158 of hMC3R, Trp9 utilizes π-π stacking interactions with F295 and F298, located on the transmembrane domain of hMC3R. It is hypothesized that as the frequency of Trp9-hMC3R interactions decrease, antagonistic activity increases. The absence of any interactions of the N-methyl groups with hMC3R suggests that their primary function is to modulate backbone conformations of the ligands.

Theoretical Evidence for Multiple Charge Transfer Pathways in Bacteriorhodopsin.

The development of molecular-scale junctions utilizing biomolecules is a challenging field that requires intimate knowledge of the relationship between molecular structure and conductance characteristics. One of the key parameters to understanding conductance efficiency is the charge mobility, which strongly influences the response time of electronic devices. The charge mobility of bacteriorhodopsin (bR), a membrane protein that has been studied experimentally in detail, was theoretically investigated using extended Marcus-Hush theory. Charge mobilities of 1.3 × 10(-2) and 9.7 × 10(-4) cm(2)/(V s) for hole and electron transfer, respectively, were determined. The computed electron mobility is comparable to experimentally measured values (9 × 10(-4) cm(2)/(V s)). Interestingly, the pathways for hole and electron hopping were very distinct from each other, utilizing different transmembrane helices to traverse the protein. In particular, only the electron transfer pathway involved the retinal chromophore, indicating that the efficiency of charge transfer is directly affected by the tertiary arrangement of proteins. Our results provide a template for obtaining the molecular and electronic-level details that can reveal fundamental insights into experimental studies on protein electron transport and inform efficient design of biomolecular-based junctions on the nanoscale.

Proteorhodopsin Activation Is Modulated by Dynamic Changes in Internal Hydration.

Proteorhodopsin, a member of the microbial rhodopsin family, is a seven-transmembrane α-helical protein that functions as a light-driven proton pump. Understanding the proton-pumping mechanism of proteorhodopsin requires intimate knowledge of the proton transfer pathway via complex hydrogen-bonding networks formed by amino acid residues and internal water molecules. Here we conducted a series of microsecond time scale molecular dynamics simulations on both the dark state and the initial photoactivated state of blue proteorhodopsin to reveal the structural basis for proton transfer with respect to protein internal hydration. A complex series of dynamic hydrogen-bonding networks involving water molecules exists, facilitated by water channels and hydration sites within proteorhodopsin. High levels of hydration were discovered at each proton transfer site-the retinal binding pocket and proton uptake and release sites-underscoring the critical participation of water molecules in the proton-pumping mechanism. Water-bridged interactions and local water channels were also observed and can potentially mediate long-distance proton transfer between each site. The most significant phenomenon is after isomerization of retinal, an increase in water flux occurs that connects the proton release group, a conserved arginine residue, and the retinal binding pocket. Our results provide a detailed description of the internal hydration of the early photointermediates in the proteorhodopsin photocycle under alkaline pH conditions. These results lay the fundamental groundwork for understanding the intimate role that hydration plays in the structure-function relationship underlying the proteorhodopsin proton-pumping mechanism, as well as providing context for the relationship of hydration in proteorhodopsin to other microbial retinal proteins.

Explaining the mobility of retinal in activated rhodopsin and opsin.

Rhodopsin, the mammalian dim light photoreceptor, is the canonical model for G protein-coupled receptors. Activation of rhodopsin occurs when the covalently bound inverse agonist, retinal, absorbs a photon and undergoes an 11-cis to all-trans isomerization. Two critical components of the visual cycle occur with the (1) hydrolytic release of all-trans retinaldehyde and subsequent (2) uptake of 11-cis retinaldehyde to reform the Schiff base linkage in the apoprotein opsin. Two pores on the surface of opsin are connected via the retinal channel, as discovered upon solution of the X-ray crystal structure (Park et al., Nature, 2008), and could serve as potential entryways for uptake and release. Using molecular dynamics simulations, we examined the behavior of rhodopsin in the Meta-II conformation (active) under Meta-I conditions (inactive), and discovered that the retinal binding pocket is flexible enough to allow a 180° rotation along the long axis of the retinal polyene chain. This result reconciles a discrepancy between the known polyene chain orientation from crystallographic and spectroscopic studies and opens the door for further investigation into the intermolecular interactions between the retinal ligand and the apoprotein opsin. Subsequent docking studies of both isomers of retinal into the opsin channel were then conducted to identify the mechanism for uptake and release. Our results suggest that retinal undergoes unidirectional uptake through Pore A and release through Pore B, and that aromatic sidechain interactions play a key role in stabilizing retinal within the opsin channel. These findings are significant in developing our understanding of the retinoid cycle and how ligand-receptor interactions in rhodopsin relate to G protein-coupled receptor activation.

Novel Phosphotidylinositol 4,5-Bisphosphate Binding Sites on Focal Adhesion Kinase.

Focal adhesion kinase (FAK) is a protein tyrosine kinase that is ubiquitously expressed, recruited to focal adhesions, and engages in a variety of cellular signaling pathways. Diverse cellular responses, such as cell migration, proliferation, and survival, are regulated by FAK. Prior to activation, FAK adopts an autoinhibited conformation in which the FERM domain binds the kinase domain, blocking access to the activation loop and substrate binding site. Activation of FAK occurs through conformational change, and acidic phospholipids such as phosphatidylinositol 4,5-bisphosphate (PIP2) are known to facilitate this process. PIP2 binding alters the autoinhibited conformation of the FERM and kinase domains and subsequently exposes the activation loop to phosphorylation. However, the detailed molecular mechanism of PIP2 binding and its role in FAK activation remain unclear. In this study, we conducted coarse-grained molecular dynamics simulations to investigate the binding of FAK to PIP2. Our simulations identified novel areas of basic residues in the kinase domain of FAK that potentially undergo transient binding to PIP2 through electrostatic attractions. Our investigation provides a molecular picture of PIP2-initiated FAK activation and introduces promising new pathways for future studies of FAK regulation.